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Polymer  electrolyte  membrane  direct  methanol  fuel  cells  (PEM-DMFCs)  have  several  advantages  over 
hydrogen-fuelled  PEM  fuel  cells;  but  sluggish  methanol  electrochemical  oxidation  and  methanol 
crossover  from  the  anode  to  the  cathode  through  the  PEM  are  two  major  problems  with  these  cells.  In  the 
present  work,  a  comprehensive  one-dimensional,  single  phase,  isothermal  mathematical  model  is  devel¬ 
oped  for  a  liquid-feed  PEM-DMFC,  taking  into  account  all  the  necessary  mass  transport  and  electrochemi¬ 
cal  phenomena.  Diffusion  and  convective  effects  are  considered  for  methanol  transport  on  the  anode  side 
and  in  the  PEM,  whereas  only  diffusional  transport  of  species  is  considered  on  the  cathode  side.  A  multi- 
step  reaction  mechanism  is  used  to  describe  the  electrochemical  oxidation  of  methanol  at  the  anode. 
Stefan-Maxwell  equations  are  used  to  describe  multi-component  diffusion  on  the  cathode  side  and  Tafel 
type  of  kinetics  is  used  to  describe  the  simultaneous  methanol  oxidation  and  oxygen  reduction  reactions 
at  the  cathode.  The  model  fully  accounts  for  the  mixed  potential  effect  caused  by  methanol  crossover  at  the 
cathode.  It  shows  excellent  agreement  with  literature  data  of  the  limiting  current  density  for  different  low 
methanol  feed  concentrations  at  different  operating  temperatures.  At  high  methanol  feed  concentrations, 
oxygen  depletion  on  the  cathode  side,  due  to  excessive  methanol  crossover,  results  in  mass-transport 
limitations.  The  model  can  be  used  to  optimize  the  geometric  and  physical  parameters  with  a  view  to 
extracting  the  highest  current  density  while  still  keeping  a  tolerably  low  methanol  crossover. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Diminishing  fossil  fuel  reserves,  ever  escalating  fuel  prices, 
contemporary  industrial  developments  and  the  rapid  pace  of 
urbanization  call  for  secured  and  environmentally  sustainable 
energy  sources.  Fuel  cells  are  currently  advocated  as  the  most 
promising  technology  to  meet  our  future  energy  requirements. 
They  offer  a  potentially  non-polluting  and  renewable  means  of 
generating  electricity.  Of  the  several  types  of  fuel  cell  available 
today,  liquid-feed  polymer  electrolyte  membrane  direct  methanol 
fuel  cells  (PEM-DMFCs)  have  been  extensively  investigated  in  the 
past  several  years  and  are  considered  as  leading  contenders  for 
transportation,  portable  and  micro-power  applications  as  they  offer 
a  combination  of  simplicity,  robustness  and  high  energy  density 
[1-3].  Nevertheless,  the  commercial  exploitation  of  liquid-feed 
PEM-DMFCs  is  being  hindered  by  two  major  technical  impedi¬ 
ments,  namely,  sluggish  electrochemical  oxidation  of  methanol 
at  the  anode,  which  results  in  high  anode  overpotentials,  and 
methanol  crossover  from  the  anode  to  the  cathode  through  the 
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PEM,  which  causes  fuel  and  oxidant  loss  (on  the  anode  and  the 
cathode  sides,  respectively)  and  cathode  electrode  depolarization. 
A  major  breakthrough  in  the  research  and  development  of  liquid- 
feed  PEM-DMFCs  has  been  the  modification  of  platinum  by  a  second 
component  such  as  Pt-Ru.  Even  with  this  binary  catalyst,  however, 
the  overpotential  at  the  anode  in  a  liquid-feed  PEM-DMFC  is  still 
much  higher  than  that  of  a  hydrogen-fuelled  PEMFCs. 

To-date,  even  for  an  elaborately  optimized  PEM  with  reasonable 
ionic  conductivity,  methanol  crossover  is  still  a  major  challenge 
[4].  Many  studies  have  been  published  in  the  literature  dealing 
with  the  issue  of  methanol  crossover,  which  is  found  to  be  strongly 
dependent  on  the  methanol  feed  concentration.  Constrained  by  the 
problem  of  methanol  crossover,  dilute  methanol  solutions  (around 
1.0  M)  are  usually  employed  in  liquid-feed  PEM-DMFCs.  Although 
further  dilution  of  methanol  solutions  will  additionally  reduce 
methanol  crossover,  mass  transport  limitation  can  become  a  seri¬ 
ous  problem  and  lead  to  lower  limiting  current  densities,  especially 
at  methanol  concentrations  less  than  0.5  M  [5]. 

Besides  the  aforementioned  phenomena,  a  variety  of  operating 
and  geometric  parameters  affect  the  performance  or  polariza¬ 
tion  behaviour  of  liquid-feed  PEM-DMFCs  [6].  Comprehensive 
experimental  investigations  on  fuel  cells  have  traditionally  been 
difficult  and  often  turn  out  to  be  practically  and  economically 
unfeasible.  Mathematical  modelling  has  proved  to  be  a  potential 
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tool  in  this  regard  to  provide  qualitative  insights  into  the  trans¬ 
port  processes  involved  in  fuel  cells.  Numerous  models  ranging 
from  one-dimensional,  single-phase  [5,7-13]  to  multiphase,  multi¬ 
dimensional  [14-17]  can  be  found  in  the  literature.  The  majority 
of  these  models  are  primarily  focused  on  overall  cell  polarization 
behaviour  without  systematically  addressing  important  phenom¬ 
ena  such  as  the  mass  transport-driven  limiting  current  behaviour, 
the  effect  of  some  key  parameters  on  methanol  crossover  and 
the  consequent  effect  on  the  cell  performance.  Further,  the  phys¬ 
ical  domain  is  simplified  by  assuming  the  catalyst  layers  to  be 
very  thin  and  the  mass-transport  limitation  on  the  cathode  side  is 
neglected.  The  electrochemical  oxidation  of  methanol  is  described 
by  Tafel  kinetics  with  a  fixed  reaction  order.  In  reality,  this  reac¬ 
tion  is  known  to  be  a  complex  multi-step  reaction  which  proceeds 
through  a  series  of  elementary  reaction  steps  [18-20].  Recent  exper¬ 
iments  by  Vidakovic  et  al.  [21]  show  that  the  reaction  order  with 
respect  to  methanol  to  varies  from  zero  (low  overpotential,  high 
concentrations)  to  close  to  unity  (high  overpotential,  low  con¬ 
centrations).  Therefore,  a  kinetic  model  must  be  employed  that 
accurately  accounts  for  the  reaction  order  transition  from  zero  to 
unity  depending  on  the  operating  conditions. 

In  the  present  work,  a  comprehensive  steady-state,  isother¬ 
mal,  one-dimensional,  single  phase,  semi-analytical  mathematical 
model  is  developed  taking  into  account  all  the  necessary  mass 
transport  and  electrochemical  phenomena  to  predict  the  perfor¬ 
mance  of  a  liquid-feed  PEM-DMFC.  The  model  also  investigates 
the  effect  of  key  operating  parameters  on  methanol  crossover  and 
methanol  mass-transport-driven  limiting  current  behaviour. 

2.  Mathematical  model 

A  liquid-feed  PEM-DMFC  operates  by  the  electrochemical  oxi¬ 
dation  of  methanol  and  water  at  the  anode  in  the  presence  of  a 
Pt-Ru/C  binary  catalyst,  i.e., 

CH3OH  +  H20  -*  C02  +  6H+  +  6e“  (E°node  =  0.02  V  vs.  SHE) 

(1) 

and  the  electrochemical  reduction  of  oxygen  at  the  cathode  in  the 
presence  of  Pt/C  catalyst: 

|02  +6H+  +6e-  -*  3H20  (Ec°athode  =  1.23  V  vs.  SHE)  (2) 

The  overall  reaction  can  therefore  be  written  as 

CH3OH  +  |02  C02  +  2H20  (E°el|  =  1.21  V)  (3) 

Fig.  1  depicts  a  schematic  diagram  of  a  liquid-feed  PEM-DMFC  seg¬ 
mented  into  seven  regions,  namely,  anode  flow  chamber  (AFC), 

CH3OH  +  HzO  +  COz  (oy  Air)  +  C02  +  HzO 


CH3OH  +  H20  (02  /  Air) 


Fig.  1.  Schematic  division  of  a  DMFC  into  seven  regions. 


Table  1 

Species  and  phenomena  considered  in  various  regions  of  DMFC. 


Region 

Species  considered 

Mass  transport 

Electrochemical 

kinetics 

AFC 

Methanol,  water,  and 
dissolved  carbon  dioxide 

V 

ABL 

Methanol,  water,  and 
dissolved  carbon  dioxide 

V 

ACL 

Methanol  and  water 

V 

^/(non-Tafel  type) 

PEM 

Methanol  and  water 

V 

CCL 

Oxygen 

V 

V  (Tafel  type) 

CDL 

Oxygen,  nitrogen,  water 
vapor,  and  carbon  dioxide 

V 

CFC 

Oxygen,  nitrogen,  water 
vapor,  and  carbon  dioxide 

V 

anode  backing  layer  (ABL),  anode  catalyst  layer  (ACL),  polymer 
electrolyte  membrane  (PEM),  cathode  catalyst  layer  (CCL),  cathode 
diffusion  layer  (CDL),  and  cathode  flow  chamber  (CFC).  The  origin 
is  set  at  the  AFC|ABL  interface.  By  convention,  any  flux  in  the  direc¬ 
tion  of  increasing  thickness  is  taken  to  be  positive.  A  summary  of 
the  species  and  the  mass  transport  and/or  electrochemical  phe¬ 
nomena  considered  in  each  of  these  regions  is  presented  in  Table  1 . 
The  model  is  constructed  based  on  the  following  assumptions. 

•  The  cell  is  operating  under  steady-state  and  isothermal  condi¬ 
tions. 

•  Variations  in  only  one  spatial  Cartesian  co-ordinate,  i.e.,  per¬ 
pendicular  to  the  membrane  electrode  assembly  (MEA),  are 
considered. 

•  The  flow  chambers  are  thoroughly  mixed. 

•  The  catalyst  layers  are  macro-homogeneous  porous  electrodes. 
Therefore,  the  electrochemical  reactions  taking  place  in  these 
regions  are  modelled  as  homogeneous  reactions. 

•  C02  is  fully  dissolved  and  no  C02  permeation  takes  place  through 
the  membrane. 

•  The  membrane  is  assumed  to  be  fully  saturated. 

•  The  electrolyte  membrane  is  impervious  to  cathode  side  gases. 

•  Only  water  vapour  is  considered  on  the  cathode  side. 

•  The  crossed-over  methanol  through  the  membrane  is  oxidized 
instantaneously  at  the  PEM  |  CCL  interface  due  to  the  large  cathode 
overpotential  imposed  on  the  methanol. 

•  A  voltage  loss  due  to  electronic  resistance  is  not  considered. 

•  The  overpotential  in  the  catalyst  layer  is  ohmic  in  nature  and 
since  the  catalyst  layer  is  very  thin,  changes  in  overpotential  are 
neglected. 

The  governing  equations,  for  the  each  of  the  aforementioned 
regions,  are  described  in  detail  in  the  following  sub-sections. 

2.1.  Mass  transport  in  the  anode  and  cathode  flow  chambers 

In  the  flow  chambers,  the  reactants  split  their  flow  between  the 
backing/diffusion  layers  and  the  channel  exits,  as  determined  by 
the  reaction  stoichiometry  and  the  respective  permeation  rates. 
Since  the  contents  in  the  flow  chambers  are  thoroughly  mixed, 
the  concentration/mole  fraction  of  the  reactants  at  the  flow- 
chamber|backing-layer  interface  is  equal  to  the  concentration/mole 
fraction  of  reactants  in  the  exit  stream,  which  can  be  determined  by 
a  mass  balance  of  the  species  across  the  flow  chambers,  as  described 
below.  For  liquid-phase  transport  in  the  AFC,  a  generalized  mass 
balance  over  the  anode-side  species  can  be  written  as 

n°f  =  n|na  -Ax  (NfBL);  i  =  Me,  W,  Q  (4) 

where  n|A ,  n°ar,  N^BL,  and  A  represent  the  inlet  and  the  outlet  molar 
flow  rates  of  species  i  on  the  anode  side,  the  flux  of  species  i  through 
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the  ABL,  and  the  active  area  of  the  cell  respectively.  Now,  since  the 
concentration  of  methanol  and  dissolved  carbon  dioxide  are  very 
low  in  the  anode-side  exit  stream,  the  mole  fraction  of  methanol  in 
the  exit  stream,  xjJJJ ,  can  be  determined  as 


^Me-w  1S  effective  diffusivity  coefficient  of  methanol  in  water 
in  the  ABL  and  is  given  by 

r»ABL,eff  /cABLd-5  r> 

UMe-\N  =  )  UMe-W 


/-OUt  /-OUt  n°Ut 

xout  _  _ ^Me _  ^  ^Me  _  Me, a  (cn 

AMe  /-out  ,  /-out  .  /-out  ~  /-out  nout  '  ' 

LW  +  LMe  +  ^Ct  LW  “a 

where  n^ra  is  the  outlet  molar  flow  rate  of  methanol  on  the  anode 
side  and  n°  ut  is  the  total  outlet  molar  flow  rate  on  anode  side,  which 
is  the  sum  of  outlet  flow  rates  of  methanol,  water  and  dissolved 
carbon  dioxide  calculated  from  Eq.  (4). 

From  the  last  equality  in  Eq.  (5),  the  concentration  of  methanol 
in  the  anode  outlet  stream,  CJJ}£a,  which  is  assumed  to  be  the  con¬ 
centration  of  methanol  at  the  AFC|ABL  interface,  CABC/ABL,  can  be 
expressed  as 


nOUt  /  _  V  ti°UL 

/-out  _  /-AFC/ABL  _  pout  nMe  _  /  A/V  \  Me, a 

Me-3  ”  “  ^W  nout  ~  \MW  J  „out 


For  gas  phase  transport  in  the  CFC,  a  generalized  mass  balance  over 
the  cathode-side  species  can  be  written  as 


=  n|nc  -  A  x  (NfDL);  i  =  0,  N,  WV,  Cg 


(7) 


where  nJA ,  n?^,  NCDL  have  the  same  meaning  as  before  except  that 
they  pertain  to  the  cathode  side.  Now,  the  mole  fraction  of  each  of 
these  species  in  the  exit  stream,  which  is  assumed  to  be  the  mole 
fraction  at  the  CFC|CDL  interface,  can  be  determined  as 


where  £ABL  is  the  porosity  of  the  ABL,  the  exponent  1.5  is  the  Brugge- 
man  correction  factor  that  accounts  for  pore  tortuosity,  and  DMe- w 
is  the  diffusivity  coefficient  of  methanol  in  water. 

Since  the  concentrations  of  methanol  and  dissolved  carbon  diox¬ 
ide  are  very  low,  Eq.  (11 )  can  be  written  as 

nrABL 

mABL  _  nABL,eff  uv~Me  ,  ,,ABL/-ABL 
iVMe  ~  ~uMe  h  vs  LMe 

where  vABL  =  (Mw/pw)NA0BL  is  the  superficial  convective  velocity 
in  the  ABL;  Mw  and  pw  are  molecular  weight  and  density  of  water, 
respectively. 

Since  there  is  no  chemical  reaction  in  this  region,  the  flux  of  any 
species  remains  constant.  Thus, 


dNABL 

—l — =  0;  i  =  Me,W,Q  (14) 

dz 

Writing  down  Eq.  (14)  for  methanol  and  combining  with  Eq.  (11), 
the  following  implicit  expression  for  methanol  flux  through  the 
ABL,  NABL,  can  be  obtained  as 


nabl  _ 

iVMe  - 


^AFC/ABL 

^Me 


expf^VfcABL  )-C^CfvABL 


exp  (vABL/kABL)  -  1 


(15) 


nout 

Yout  _  yCFC/CDL  _ 

Xi,c  ~Xi  ~  ^out  W 

where  n?^  is  the  outlet  molar  flow  rate  of  species  i  on  cathode  side, 
and  n°ut  is  the  total  molar  flow  rate,  which  is  the  sum  of  molar 
flow  rates  of  oxygen,  nitrogen,  water  vapour  and  carbon  dioxide 
calculated  from  Eq.  (7)  on  the  cathode-side,  respectively. 


where  /<ABL  =  DABL,eff/^ABL  is  the  mass  transfer  coefficient  of 
methanol  in  the  ABL,  8ABl  is  the  thickness  of  the  ABL,  and  CABL/ACL 
is  the  concentration  of  methanol  at  the  ABL|ACL  interface. 

2.3.  Mass  transport  and  electrochemical  phenomena  in  the  anode 
catalyst  layer 


2.2.  Mass  transport  in  anode  backing  layer 


A  mixture  of  methanol,  water  and  dissolved  carbon  dioxide  are 
considered  in  this  region.  Dissolved  carbon  dioxide  diffuses  back 
into  the  AFC  and  exits  with  the  solution.  The  flux  of  each  of  these 
species  is  composed  of  two  factors:  one  required  for  the  reaction 
and  one  that  diffuses  (crosses  over)  through  the  membrane.  There¬ 
fore  the  flux  in  this  region  is  expressed  as 

NfBL  =  Nfto  +  NfEM ;  i  =  Me,  W,  C,  (9) 


where  N?t0  and  /Vr’EM  represent  the  stoichiometric  flux,  i.e.,  flux 
required  for  the  reaction  and  the  flux  diffusing  through  the  PEM 
(crossover  flux),  respectively,  of  species  i. 

The  stoichiometric  flux  is  determined  by  the  reaction  stoichiom¬ 
etry: 


ivrSto 

iVMe 


ivrStO  _  _  luStO 
iVW  - 


I 

6 F 


(10) 


where  /  is  the  cell  current  density  and  F  is  the  Faraday  constant 
(96,  485  C  mol-1).  As  assumed  earlier,  the  crossover  flux  of  dis¬ 
solved  carbon  dioxide  is  set  to  zero.  Methanol  transport  in  this 
region  is  a  combination  of  diffusion  and  convection.  Therefore,  the 
flux  of  methanol,  NABL,  is  expressed  as 


Ni 


rABL 


nABL,eff 

'^Me-W 


dCABL 

ULMe 

dz 


■  vABLivrABL 
+  AMe  iVTot 


(ii) 


where  x^BL,  NABtL  are  the  local  mole  fraction  of  methanol  in  the  ABL, 
and  total  flux  through  the  ABL,  respectively.  Here 


<BtL  =  NMeL  +  NWL  +  NC1BL 


(12) 


2.3 A.  Mass  transport  in  anode  catalyst  layer 

A  fraction  of  the  methanol  solution,  arrived  in  this  region 
through  the  ABL,  gets  electrochemically  oxidized  in  the  presence 
of  Pt-Ru/C  catalyst  via  the  electrochemical  reaction  shown  in  Eq. 
(1 )  and  the  rest  is  transported  through  the  PEM  by  different  driving 
forces  depending  upon  the  operating  conditions. 

A  binary  mixture  of  methanol  and  water  is  considered  in  this 
region  and  the  flux  of  each  of  these  species  is  assumed  to  vary 
linearly  as 


ABL 


JUALL  _  jui 

iVMe  “iVMe 


/ 

6F 


z  _ZABL/ACL 
SACL 


(16) 


J\jALL  _  iWi 

1  ’  \  A/  —  1  n 


w 


/ 

6F 


z  _ZABL/ACL 
SACL 


(17) 


where  8ACl  is  the  thickness  of  the  ACL. 

Therefore,  the  total  flux  variation  in  the  ACL,  NA^tL,  is  given  by 


mACL  _  iuABL  _ 
iVTot  “  iVTot 


/ 

3 if 


z  _ZABL/ACL 
8 


(18) 


Methanol  transport  in  this  layer  is  a  combination  of  diffusion  and 
convection.  Therefore,  the  flux  of  methanol  in  this  region  N^L  is 
expressed  as 


Ni 


Me 


nACL,eff 

'^Me-W 


dcACL 

UL-Me 

dz 


+  NTotL*MeL 


(19) 


where  NA^tL  is  the  averaged  total  flux  obtained  from  Eq.  (17)  and 
^Me-w  1S  effective  diffusivity  coefficient  of  methanol  in  water 
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in  the  ACL  and  is  given  by 


For  the  reaction  step  (22c),  the  rate  expression  is  written  as 


nACL,eff 

^Me-W 


(£ACL)  15DMe-W 


where  £ACL  is  the  porosity  of  the  ACL. 
Eq.  (19)  can  be  further  written  as 


N; 


Me 


nACL,eff 

'^Me-W 


dCACL 

ULMe 

dz 


(20) 


where  t£CL  =  (Mw/pwJM^-1  is  the  average  superficial  convective 
velocity  in  the  ACL. 

Equating  Eqs.  (16)  and  (20)  and  upon  solving  with  the  boundary 
condition  CA^L  =  CA®L/ACL  at  z  =  zABL/ACL,  the  following  expression  for 
the  distribution  of  methanol  concentration  in  the  ACL,  C^L(z),  can 
be  obtained,  i.e., 


rACLr~>i  _  rABL/ACL  f  ^CL  ^ABL/ACL 

^Me  W  ^Me  ^  l  nACL,eff  ^Me 

V^Me-W 


NABL  , 

Me  \  (z  _ZABL/ACL) 


,ACL,eff 

Me-W 


(21) 


2.3.2.  Kinetic  model  for  methanol  electrochemical  oxidation  in 
anode  catalyst  layer 

A  simplified  Gasteiger  mechanism  [22]  is  used  in  the  present 
work  to  develop  a  kinetic  expression  similar  to  that  of  Meyers  and 
Newman  [23]  describing  the  reaction  rate  of  methanol  in  terms  of 
methanol  concentration  and  the  anode  overpotential,  and  is  later  on 
used  to  obtain  anode  overpotential  data.  The  reaction  mechanism 


consists  of  the  following  elementary  steps: 

CH3OHafc  -»  CH3OHACl  (mass  transfer)  (22a) 

CH3OH  +  Pt  ->  CH3OHadPt  ->  •  •  COadPt  +  4H+  +  4e“  (22b) 

H20  +  Ru  4^  OHadRu  +  H+  +  e-  (22c) 

COadPt  +  OHad  Ru  — C02  +  H+  +  e  +  Pt  +  Ru  (22d) 


As  a  first  step  (22a),  methanol  is  transferred  from  the  AFC,  through 
the  ABL,  into  the  ACL  wherein  the  electrochemical  oxidation 
of  methanol  takes  place  through  several  steps.  The  potential- 
dependent  final  step  (22d)  involving  the  reaction  between  two 
surface-adsorbed  species  is  the  rate-determining  step,  except 
at  sufficiently  high  current  densities,  where  the  first  step 
(potential-independent  methanol  physiosorption)  of  the  consec¬ 
utive  irreversible  reaction  sequence  (22b)  limits  the  reaction.  The 
third  step  (22c),  i.e.,  the  water  dissociative  adsorption,  is  assumed 
to  be  in  equilibrium. 

The  balance  of  adsorbed  species  on  the  surface  of  the  two  cata¬ 
lyst  components,  platinum  and  ruthenium,  yields: 


$pt  +  $c,Pt  =  1  (23  a) 

where  0c,pt  refers  to  fractional  coverage  on  Pt  by  some  C-containing 
adsorbates  without  specifying  the  type  of  the  adsorbed  species. 
Similarly: 

Oru  +  #OH,Ru  =  1  (23b) 

The  rate  expressions  for  every  step  considered  in  the  above  mecha¬ 
nism  are  formulated  as  follows.  The  rate  expression  for  the  reaction 
step  (22b)  is  written  as 

ri=/<iCMe0pt  (24) 

where  r\  and  k\  represent  the  rate  of  reaction  and  the  reaction  rate 
constant  of  step  (22b),  respectively. 

From  Eqs.  (23a)  and  (24),  we  have 

ri=/<iCMe(l-0c,Pt)  (25) 


r2  =  /<2)faw0Ru  exp 


a2FV a- 
RT 


/<2,b^OH,Ru  exp 


r-(l  -a2)FVal 

RT 


(26) 


where  r2,  k2 ,f,  k2ib,  and  ol2  represent  the  rate  of  reaction,  the  forward 
and  the  backward  chemical  reaction  rate  constants,  and  the  transfer 
coefficient  of  step  (22c),  respectively,  and  Va  is  the  anode  voltage. 
From  Eqs.  (23b)  and  (26),  we  have 


*2  =  fc2,a0w(l-#OH,Ru)eXp 


'a2FVa' 

-^2,b^OH,Ru  exp 

r— (i— «2)FVal 

RT  ] 

RT 

(27) 

where  aw  is  the  activity  of  water. 

As  this  reaction  is  assumed  to  be  in  equilibrium,  Eq.  (27)  can  be 
solved  for  0Oh,ru  as  (water  activity,  aw  is  unity): 


a  I<2  exp (FVa/RT) 

0H'Ru  1  +  l<2  exp  (FVa/RT) 


(28) 


where  I<2  =  k2  flk2  b  is  the  chemical  equilibrium  constant  of  reaction 
step  (22c). 

For  the  reaction  step  (22d),  the  rate  expression  is  written  as 


r3  =  k30c,  Pt^oH.Ru  exp 


'a  3FVa' 
RT 


(29) 


where  r3,  k3  and  a3  represent  the  rate  of  reaction,  the  chemical 
reaction  rate  constant  and  the  transfer  coefficient  of  step  (22d), 
respectively. 

Applying  steady-state  conditions,  i.e.,  r\  =r3,  one  can  obtain  an 
expression  for  0C, pt  as 

a  _  k1CMe(l+^2exp(FVa/RT)) 

C'Pt  k3l<2  exp((a3  +  1  )FVa/RT)  +  fc,  CMe(l  + 1<2  exp (FVa/RT)) 

(30) 


Now,  the  overall  reaction  rate  can  be  expressed  by  either  the  rate 
of  reaction  of  step  (22b),  or  by  the  rate  of  reaction  of  step  (22d). 
Since  six  electrons  are  exchanged  in  the  oxidation  reaction,  the 
corresponding  anodic  current  density,  ja,  can  be  obtained  as 


j  a  =  6Fr3  =  6F/<30c,Pt0OH,Ru  exp 


cx3FVa 

RT 


(31) 


Combining  Eqs.  (27),  (29)  and  (30)  and  assuming  I<2  exp (FV/RT) » 1, 
we  get  for  the  anodic  current  density  the  following  expression: 

•  _  fiF  fe3CMe  exp(g3FVa/RT) 

Cue  +  fe/ki )  exp(a3FVa/Rr) 

Replacing  a3  by  aa,  the  anodic  current  density  given  in  Eq.  (32)  can 
be  further  expressed  in  terms  of  anode  overpotential  as 

•  =  PCue  exp(aai]aF/RT) 

CMe  +  Y  exp(aaF?)a/RT) 

The  kinetic  parameters  appearing  in  Eq.  (33),  namely,  fi  and  y,  are 
given  by 


/?  =  6Fk3exp(^) 

k3  faaFEa\ 

y=k-,exp{-TT  ) 


where  Fa  is  the  equilibrium  anode  potential. 

The  transition  of  the  reaction  order  from  zero,  at  low  overpo¬ 
tentials  (the  second  term  in  the  denominator  of  Eq.  (33)  can  be 
ignored,  i.e.,  the  reaction  is  potential  dependent),  to  one  at  high 
overpotentials  (concentration  of  methanol  becomes  low;  therefore 
the  first  term  in  the  denominator  of  Eq.  (33)  can  be  ignored,  i.e., 
the  reaction  is  concentration  dependent)  can  be  observed  from  Eq. 
(33). 
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Integration  of  Eq.  (33)  over  the  anode  catalyst  layer  yields  the 
cell  current  density,  /,  as 


current  density,  J,  and  pressure  gradient  across  the  PEM,  dp/dz,  as 
shown  in  Eq.  (39): 


j  =  I**1  /^(Z)exp(«a>?aF/fir)  dz 
Jo  *-MeL(z)  +  y  sxp(aaFr/a/RT) 

This  kinetic  expression  along  with  Eq.  (21)  is  used  to  calculate  the 
anode  overpotential  for  any  given  cell  current  density. 

2.4.  Mass  transport  in  polymer  electrolyte  membrane 


d«/>  0/F)  +  C^+M(fcp/ /x)(dP/dz) 

dz  ~  FCPEMaDPEM/RT)  +  (k0/^)CPEM) 

Substituting  Eq.  (39)  in  Eq.  (36),  an  explicit  expression  for  the  pore 
convective  velocity  in  the  PEM,  vPEM,  can  be  derived  as 


k0  r(//F)  +  CPEM(kp//r)(dp/dz)" 
77  |_((DPEM /RT)  +  (k<p//r)CPEM)_ 


kp  dp 
/x  dz 


(40) 


Following  the  approach  of  Bernardi  and  Verbrugge  [24],  the 
Nernst-Planck  equation  with  the  modified  Schlogl  equation  is  used 
to  characterize  the  transport  of  species  in  the  membrane  electrolyte 
region: 


Water  transport  across  the  PEM  takes  place  via  diffusion  and  con¬ 
vection.  Since  the  membrane  is  assumed  to  be  fully  saturated,  only 
convective  water  transport  is  considered.  Therefore,  the  flux  of 
water  across  the  PEM,  N^M,  is  given  by 


N PEM  _  2pem  r  £)PEM^PEM^0  _  £)PEM  +  ^PEM  yPEM 

~~  4  1  dz 


F 

Rf^ 


dz 


(35) 


In  accordance  with  the  modified  Schlogl  equation,  the  pore  convec¬ 
tive  velocity,  vPEM,  is  given  by 


l/PEM  =  !^Z1 


/-lJt.lVl  p  4P 

f  dz  fi  dz 


7PEMrPEMi:^ 


(36) 


where  k0  is  the  electrokinetic  permeability,  CPEM  is  the  concentra¬ 
tion  of  fixed  charge  species  (sulfonic  acid  groups)  in  the  PEM,  ZPEM 
is  the  charge  number  of  fixed  charge  species  in  the  PEM,  kp  is  the 
hydraulic  permeability,  and  pi  is  the  liquid  water  viscosity. 

For  neutral  species,  the  first  term  of  Eq.  (35)  vanishes.  Therefore, 
the  superficial  flux  of  neutral  species  in  the  membrane  electrolyte 
region  is  expressed  as 


d^PEM 

Npem  =  -DPEfM  — \ —  +  cPEMvPEM;  i  =  Me,  W 

1  1,011  /-J7  l 


(37) 


mPEM  _  rPEM.,PEM 
iVMe  -  LW  ys 


VPEM 


Vw 


(41) 


Methanol  is  transported  through  the  PEM  in  the  same  way  as  it 
is  transported  through  the  ABL.  So,  the  flux  of  methanol  through 
the  PEM  or  the  crossover  flux  of  methanol  can  be  expressed  in  a 
similar  form,  as  shown  in  Eq.  (15).  However,  since  the  crossed-over 
methanol  is  assumed  to  be  instantaneously  oxidized  at  the  PEM|CCL 
interface,  this  equation  can  be  further  reduced  to  the  form  shown  in 
Eq.  (42)  by  setting  the  concentration  of  methanol  at  the  PEM/CCL, 
C-/pEM  equal  to  zero: 


npem 

iVMe 


C/PEMex  p^M/fcPEM) 
exp(vPEM/kPEM)-l  s 


(42) 


Six  protons  and  six  electrons  are  released  for  each  mole  of  methanol 
consumed  by  the  cathodic  methanol  oxidation  reaction.  Therefore 
the  crossover  current  density,  JCross,  can  be  calculated  by  multiply¬ 
ing  Nj™  by  a  factor  of  6F,  as  shown  in  Eq.  (43): 


where  vPEM  is  the  superficial  convective  velocity  in  the  membrane 
and  is  given  by 

..PEM  _  CPEM.,PEM 
vs  ~  fcW  v 

Here  £^EM  is  the  volume  fraction  of  water  in  the  PEM  which  is  given 
by 

gPEM  _  _ ^ _ 

w  ( yPEM/yw)  +  x 

where  vPEM,  Vw  and  A  are  the  partial  molar  volume  of  the  mem¬ 
brane,  the  partial  molar  volume  of  water  and  the  membrane  water 
loading  or  water  content,  respectively. 

The  equivalent  weight  (EW)  of  the  membrane  (grams  of  dry 
polymer/moles  of  sulfonic  acid  groups)  and  VPEM  are  interrelated: 

yPEM  ^ 

Po 

where  po  is  the  density  of  the  dry  membrane.  For  Nafion®  117, 
115,  or  112,  the  membranes  of  interest  here,  EW=1100,  and 
Po  =  2050kgm-3.  The  water  uptake  values  of  Nafion  from  liquid 
water  determined  by  Hinatsu  et  al.  [25]  are  used  in  the  present 
work. 

As  no  chemical  reaction  takes  place,  the  flux  of  methanol  and 
water  remain  constant: 

dNPEM 

— f —  =  0;  i  =  Me,  W  (38) 

dz 

Inserting  Eq.  (36)  into  Eq.  (35)  and  by  applying  it  for  proton  trans¬ 
port,  an  expression  for  the  potential  gradient  across  the  PEM,  d0/dz, 
can  be  formulated  in  terms  of  two  known  parameters,  namely,  cell 


f  cross  =  6FNj^EeM  (43) 

Once  the  crossover  current  density  is  calculated,  the  methanol 
crossover  efficiency  and  the  Faradaic  or  fuel  utilization  efficiency, 
JCross  and  /Eardaic,  respectively,  can  be  calculated  in  percentage  terms 
as  follows: 

rCross 

fross  {%)=wo__  (44) 

/ Fardaic  (%)=100_^_  (45) 

2.5.  Mass  transport  in  cathode  diffusion  layer 


The  gas  phase  in  the  CDL  of  a  liquid-feed  PEM-DMFC  is  a  mix¬ 
ture  of  oxygen,  nitrogen  (if  air  is  used  as  cathode  feed),  carbon 
dioxide,  and  water  vapour.  The  gaseous  mixture  is  assumed  to  be 
an  ideal  gas.  Only  diffusional  effects  are  considered  in  this  region 
and  the  Stefan-Maxwell  equations  are  used  to  account  for  multi- 
component  diffusion.  For  the  diffusion  of  an  n-c omponent  ideal  gas 
through  a  porous  medium,  the  one-dimensional  Stefan-Maxwell 
equations  take  the  form: 

HyCDL  n 

i  =  1, 2, . . . ,  n  (46) 

dz  ^  cDe  ff  J  J  1  ' 

j= i  H 


where  cDf^.  is  the  effective  concentration  binary  diffusion  coeffi¬ 
cient  and  is  given  by 


cDL 


,CDL,eff 


i-J 


=  (£CDL)  CD, 


l-J 


where  cDt_j  is  the  concentration  binary  diffusivity  coefficient  and 
£cdl  is  the  porosity  of  the  CBL. 
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Nitrogen  does  not  participate  in  any  of  the  chemical  reactions; 
it  is  inert.  Therefore,  under  steady-state  conditions,  it  remains 
motionless  and  hence  its  flux  is  set  to  be  zero: 


n£dl  =  0 


(47) 


The  fluxes  of  the  other  three  species  are  obtained  as  described 
below. 

The  flux  of  oxygen  in  this  region  is  made  up  of  two  factors:  the 
flux  required  for  the  cathode  half-cell  reaction  and  that  required  for 
the  complete  oxidation  of  the  crossed-over  methanol.  Therefore, 
the  flux  of  oxygen  through  the  CDL,  NqDL,  and  is  given  by 


(-L 

\4F 


=  -l—+-N 


vMe 


(48) 


Only  water  in  vapour  form  is  considered  in  this  region.  The  flux  of 
water  vapour  is  made  up  of  three  factors:  the  flux  produced  by  the 
cathode  half-cell  reaction,  that  produced  by  the  complete  oxidation 
of  the  crossed-over  methanol,  and  the  flux  of  water  through  the 
PEM,  N™1.  Therefore,  the  flux  of  water  vapour  through  the  CDL, 
N^L,  is  given  by 

=  i  +  2NMeM  +  N^EM  (49) 

One  mole  of  carbon  dioxide  will  be  formed  for  every  mole  of 
crossed-over  methanol  consumed.  Therefore,  the  flux  of  carbon 
dioxide  through  the  CDL,  NjJL,  is  given  by 


nP-  =  n™ 


(50) 


Since  no  chemical  reaction  occurs  in  this  region,  the  flux  of  each  of 
these  species  remains  constant  leading  to 


dNCDL 

dz 


=  0;  i  =  0,  WV,  Cg 


Applying  Eq.  (46)  for  nitrogen  (N),  oxygen  (0)  and  water  vapour 
(WV),  the  following  set  of  linear  differential  equations  can  be 
obtained: 


dx^DL 

dz 


NCDL  .wCDL  N^DL 

iVo  ,  iVwv  ,  cg 


cDeff  cDeff  cDeff 

tUN-0  tUN-W  tUN-C 


VCDL 


(52a) 


*r(z)  =  j-(mi  -41)exp(m1z)+ ^-(m2 -7ii)exp(m2z) 
B\ 


£l(Kl  —  A] )  —  C, 


B  i 


XmBL(z)  + 


B,cDf_c 


A1E2 

B , 


(53b) 


where 

x^(z)  =  x“L/CFCexp(K1(z-zCDL/CFC)) 


(53c) 

The  constants  appearing  in  the  above  equations  are  given  below: 


/  NnDl 
—  I ff — f 


^2 


B  i 


cDeff 

LUo_c 


(  n^l 

Wwv-c 

=  (a r 


+ 


\  cneff 
V  LUo-c 


B2  = 


Ng°L 

,  C^WV-0 


MCDL 

Jvwv 

cDeff 

LUo-wv 

MCDL 

Jvwv 

cDeff 

LUWV-0 


JVg°L 

cDeff 
LUo-wv 

MCDL 

,  wv 
cDeff 

LUWV-C 


c^wv-c  , 


Ci  = 


( NoDl  NoDl 


cDeff 

tU0-N 


C2  = 


MCDL 

nwv 


MCDL 

iVwv 


cDeff  cDeff 

LUWV-C  tUWV-N 


/Ci  = 


n£dl 

cDeff 

tUN-0 


ncdl  N^dl 

nWV  +  cg 


cDeff 

LiyN-WV 


cDeff 

CiyN-C 


(51)  Ei  = 


£2  = 


(C1/C1  -B2Ci  +  Ei  C2) 

*?-(*!  +B2)K]  +(A]B2-A2B]) 

(N^DLB2/cDf_c)  -  (N“lB,  /cD^_c) 


(mi,m2)  = 


AiE2-^2^i 

(Ai  +  B2)  ±  \J (Ai  +  B2 )2  -  4(Ai B2  -  A2B\ ) 


dx£DL 

dz 


NSdl  NfDL 


Lg 


■jeff  + 

'o-C 

cDeff  + 

LU0-W 

cDeff 

LUo-c 

JVg°L 

NgDL  ' 

\  yCDL 

1  xwv 

cDeff 

LUo-c 

cDeff 
LUo-w . 

ncDL  > 

d*wv 

dz 


V^S- 


nSdl 


0— N  / 


MCDL 

Jvwv 


Ng°L 

cDeff 

LUo-c 


NPl 


cd: 


.eff 


cd: 


,eff 


cd: 


,eff 


+ 


+ 


MCDL 

Jvwv 

cDeff 

LUWV-C 

MCDL 

Jvwv 

cDeff 

LUwv-c 


MCDL 

Jvwv 

cDeff 

Li^wv-0 

MCDL 

Jvwv 

cDeff 

L1A/W-N 


VCDL 


yCDL 

*wv 


MCDL 

iVwv 

cDeff 

Li7wv-c 


(52b) 


(52c) 


Eqs.  (52a)-(52c)  can  be  solved  analytically  for  the  species  concen¬ 
tration  distributions  across  the  CDL,  and  the  resulting  equations  are 
given  in  Eqs.  (53a)-(53c). 


VCDL, 


(z)  =  Ci  exp(miz)  +  c2  exp(m2z)  +  EiX^UL(z)  +  E2 


(53a) 


ci  = 


c2  = 


(m2-A-[  )XqDL/cfc 
-E2m2+{N^L/cDf c) 


-(E1m2-£1/C1  +  C,)x“ 


.CDL/CFC 


(m2  -  mi)exp(miZCDL/CFC) 


(m,  -A,)x§dl/c fc 
-£2mi  +(N“L/cD^cj 


(mi  -  m2)exp(m2zCDL/CFC) 


The  binary  concentration  diffusivity  coefficients,  cDt_j  are  esti¬ 
mated  using  the  following  correlation  [26]: 


cDab  =  2.2646  x  10 


,-5 


Y— +  — V 


1 


\Ma  Mb  7  (j2 Rlf2D)AB 


The  dimensionless  quantity  ^2D  AB  -  the  “Collisional  integral”  for 
diffusion  -  is  a  function  of  the  dimensionless  temperature  /<T/cAB. 
The  collision  integrals  have  been  curve  fitted  by  Neufeld  et  al.  [27] 
as  follows: 


_  1.06036 

“^D,AB  T,.n  iccm  + 


0.19300 


1.03587 


T*o.i56io  exp(0.47635T*)  exp(1.52996T*) 

1.76474 


exp(3. 89411 T* 
where  7^  =  /<T/cAB. 
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Here,  £AB  =  v^a^b  and  <rAB  =  (1/2 )(<ta  +  crB)  are  the  molecu¬ 
lar  separation  at  collision  and  the  energy  of  molecular  attraction 
respectively  [28]. 

2.6.  Mass  transport  and  electrochemical  phenomena  in  cathode 
catalyst  layer 

2.6 A.  Mass  transport  in  cathode  catalyst  layer 

Oxygen  transport  in  this  region  is  governed  by  diffusion.  There¬ 
fore  the  flux  of  oxygen  in  this  region,  NqCL,  is  given  by  Fields  law 
as 


drCCL 

u  °  dz 


(54) 


where  is  the  effective  diffusion  coefficient  of  oxygen  in  the 

CCL  and  is  given  by 


,CCL 


D£CLeff  =  (eCCL)  '  Dq 

where  eCCL  is  the  porosity  of  the  CCL  and  DqCL  is  the  diffusivity 
coefficient  of  oxygen  in  the  CCL. 

Applying  material  balances  based  on  standard  porous-electrode 
theory,  we  have 


d  N£CL 
dz 


so  dje 
ncF  dz 


(55) 


where  jc,  s o  and  nc  are  the  cathodic  current  density,  the  stoichio¬ 
metric  coefficient  of  oxygen  and  the  number  of  electrons  exchanged 
in  the  cathode  half-cell  reaction,  respectively. 

2.6.2.  Kinetic  model  for  electrochemical  oxygen  reduction  in 
cathode  catalyst  layer 

The  electrochemical  reduction  of  oxygen  in  this  region  is 
described  by  Tafel  kinetics  with  first-order  dependence.  Therefore, 
the  cathode  current  density  is  given  by 


Ok  -c/o 
dz  _  °'ref 


CSCL(z) 


Co, ref 


exp 


(acric{z)F^ 


(56) 


where  a/°  ref  is  the  cathode  volumetric  reference  exchange-current 
density,  pc  is  the  cathode  overpotential,  and  ac  is  the  cathode  trans¬ 
fer  coefficient. 

From  Eqs.  (55)  and  (56),  the  following  differential  equation  can 
be  derived: 

d2C^(z) 


dz2 

where 

§2  = 


—  ?  Q)  (z)  =  o 


(57) 


Soalg 


Lo^ef 


p- - exp(acr/cF/RT) 


Solving  Eq.  (57)  with  the  following  boundary  conditions: 
CSCL  =  C™L  at  z  =  zCCL/CBL 


dcCCL 

0  at  z  =  zPEM/CCL 
dz 


yields  an  expression  for  the  oxygen  distribution  in  the  CCL,  namely: 

(58) 


cccl(z)  =  CqCl/cbl  cosh^z  ~  zPEM/CCL)] 


cosh[^ca] 

Integration  of  Eq.  (56)  over  the  cathode  catalyst  layer  gives  the  cell 
current  density,  /,  as 

^CCL,eff  £2  rSCCL 


jeff  = 


ncFD ^ 


So 


/ 


Cq(z)6z 


(59) 


where  /eff  =  /+ jCross.  Thus,  the  mixed  potential  effect  on  the  cathode 
side  is  taken  into  account  by  replacing  /  by  /eff,  which  is  the  sum  of 
the  electrical  current  density  and  the  crossed-over  current  density, 
JCross,  attributed  to  the  oxidation  of  crossed-over  methanol  in  the 
CCL.  JCross  can  be  obtained  from  Eq.  (43)  in  Section  2.4.  As  will  be 
shown  later,  this  approach  correctly  captures  the  reduction  in  the 
open-circuit  voltage  at  high  methanol  concentrations. 

Solving  Eqs.  (58)  and  (59)  together,  Jeff,  takes  the  following  form: 


jeff  = 


ncfDg^effg/ 
s0 


.CCL/CBL 

^0 


tanh(^CCL) 


(60) 


Eq.  (60)  can  be  used  to  calculate  the  cathode  overpotential  for  any 
given  cell  current  density. 

2.7.  Overall  cell  voltage 

The  overall  cell  voltage  is  calculated  employing  the  following 
expression: 

£PEM 


V  =  Ece\\-rh-hc-I- 


rPEM 


(61) 


where  ECeli  is  the  difference  between  the  half-cell  potentials  of  the 
anode  and  the  cathode  corrected  for  the  thermodynamic  effects  of 
temperature  (T)  and  pressure  (P)  as  follows: 


ECen=E°+(T-T0) 


ANRT 


9£ceiA  _ 
dT  J  ncF 


In 


(k) 


(62) 


where  £°ell  is  the  ideal  electromotive  force  under  standard  condi¬ 
tions  of  temperature  (T0)  and  pressure  (P0).  (3ECell/3T)  is  the  rate 
of  change  of  ECell  with  T  and  the  last  term  represents  the  effect  of 
pressure  on  cathode  potential.  Since  water  vapour  is  assumed  to 
be  produced,  AN  is  taken  to  be  +0.5.  The  third  term  in  Eq.  (61) 
represents  the  overpotential  due  to  ohmic  resistance. 

3.  Results  and  discussion 

The  above  semi-analytical  mathematical  model,  consisting  of 
several  non-linear,  coupled  algebraic  equations,  has  been  solved 
under  galvanostatic  conditions  (cell  current  density  is  fixed  to 
obtain  the  corresponding  cell  voltage)  using  MATLAB®  over  the 
range  of  physiochemical,  geometric  and  kinetic  parameters  listed 
in  Tables  2a  and  2b.  Most  of  these  parameters  are  taken  from  the 
literature,  as  indicated.  Transport  parameters  such  as  the  diffusiv- 


ity  coefficient  of  methanol  in  water  on  the  anode  side  and  the 
binary  concentration  diffusivity  coefficients  on  the  cathode  side 

Table  2a 

Geometric,  physicochemical  and  operating  parameters  used  in  the  calculations. 

Active  area  of  the  cell  (A) 

9  cm2  [30] 

Thickness  of  the  anode  backing  layer  (5ABL) 

0.03  cm  [30] 

Thickness  of  the  anode  catalyst  layer  (<5ACL) 

0.001  cm  (assumed) 

Thickness  of  the  membrane  (<5PEM) 

0.01778  cm  (Nation  117) 

Porosity  of  the  anode  backing  layer  (£ABL) 

0.65  (assumed) 

Porosity  of  the  anode  catalyst  layer  (£ACL) 

0.4  (assumed) 

Thickness  of  the  cathode  backing  layer  (<5CBL) 

0.03  cm  [30] 

Thickness  of  the  cathode  catalyst  layer  (<5CCL) 

0.001  cm  (assumed) 

Porosity  of  the  cathode  backing  layer  (£CBL) 

0.4  (assumed) 

Porosity  of  the  cathode  catalyst  layer  (£CCL) 

0.3  (assumed) 

Air  inlet  flow  rate 

600  ml  min-1  (assumed) 

Electrokinetic  permeability  ( k 0) 

1.13  x  10-15  cm2  [24] 

Hydraulic  permeability  (fcp) 

1.58  x  10-14  cm2  [24] 

Anode  kinetic  parameter  (ft) 

6.45  x  10-6  (calibration) 

Anode  kinetic  parameter  (y) 

3.086  x  10~7  (calibration) 

Anode  transfer  coefficient  (0(3) 

0.45  (calibration) 

Cathode  transfer  coefficient  (ac) 

0.8  [14] 

Reference  cathode  exchange  current  density  ( JB  f ) 

3.75  x  10-7  Acm-2 

(calibration) 

Concentration  of  protons  in  the  membrane  (CPEM) 

1.2  x  10~3  mol  cm-3  [24] 
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Table  2b 

Temperature  dependant  physicochemical  parameters. 


Property  (K) 

343 

348 

353 

358 

363 

Diffusivity  coefficient  of  methanol  in  water  DMe_w  (cm2  s_1 ) 

3.3  x  10-5 

4.5  x  10~5 

4.1  x  10-5 

4.5  x  10-5 

5.3  x  10-5 

Diffusivity  coefficient  of  methanol  in  the  membrane  (cm2  s-1 ) 

6.06  x  10-6 

6.72  x  10~6 

7.42  x  10~6 

8.17  x  10~6 

ID 

1 

o 

X 

o 

cri 

Membrane  water  content  A  [25] 

15.4 

16.0 

17.0 

17.5 

20.0 

Viscosity  of  liquid  water,  fi  (g  cm-1  s-1 ) 

4.061  x  10-3 

3.8  x  10-3 

3.565  x  10-3 

3.55  x  10-3 

3.165  x  10“3 

are  estimated  making  use  of  empirical  correlations  available  in  the 
literature.  Results  from  the  calculations  are  discussed  below. 

3.1  Limiting  current  for  low  methanol  feed  concentrations 

For  liquid-feed  PEM-DMFCs  operating  with  low  methanol  feed 
concentrations,  limiting  current  behaviour  can  be  observed  at  rel¬ 
atively  low  current  densities.  This  behavior  can  be  attributed  to 
mass-transport  limitations  of  methanol  supply  to  the  reactive  zones 
in  the  anode  compartment  [29].  Therefore,  the  limiting  current  den¬ 
sity  may  be  interpreted  as  the  current  density  that  corresponds  to 
zero  methanol  concentration  in  the  ACL  or  as  the  current  density 
at  which  the  flux  of  methanol  through  the  ABL  equals  the  rate  of 
electrochemical  oxidation  of  methanol  in  the  ACL.  This  can  be  seen 
in  Fig.  2  where  the  predicted  variation  of  methanol  concentration  in 
the  ABL,  ACL  and  PEM  is  shown  for  different  current  densities  and 


Fig.  2.  Predicted  methanol  concentration  profiles  in  the  anode  compartment  and 
the  PEM  regions  for  a  methanol  feed  concentration  of  (a)  0.125  M,  (b)  0.25  M  and 
(c)  0.5  M  at  a  methanol  feed  flow  rate  of  1.36  ml  min-1  at  90  °C,  and  at  an  air  feed 
cathode  pressure  of  2  bar. 


for  different  methanol  feed  concentrations.  In  Fig.  2a,  for  a  methanol 
feed  concentration  of  0.125  M  at  90  °C  and  at  a  methanol  feed 
flow  rate  of  1.36  ml  min-1,  the  methanol  concentration  decreases 
in  the  ACL  as  the  current  density  increases  and  becomes  zero  at 
7  =  0.045  A  cm-2.  Decreasing  methanol  concentration  with  increas¬ 
ing  current  density  can  also  be  observed  at  the  AFC|ABL  interface. 
For  feed  concentrations  0.125  and  0.5  M,  and  at  the  same  tempera¬ 
ture  and  methanol  feed  flow  rates,  this  happens  at  current  densities 
of  0.093  and  0.197  A  cm-2,  respectively  (see  Fig.  2b  and  c).  It  can 
be  seen  from  these  methanol  concentration  profiles  that  a  higher 
flux  of  methanol  in  the  anode  catalyst  layer  cannot  be  sustained. 
This  constitutes  mass-transport  limitation  on  the  anode  side  which 
limits  the  current  density  that  can  be  extracted. 

The  limiting  current  densities  predicted  by  the  present  model 
for  three  different  methanol  feed  concentrations  (0.125,  0.25  and 
0.5  M)  at  70  and  90  °C  are  compared  with  experimental  data  taken 
from  literature  [10,30]  in  Fig.  3a  and  b,  respectively.  Further  com¬ 
parison  of  the  predicted  and  the  experimental  limiting  current 
densities  [30]  is  shown  in  Fig.  3c  as  a  function  of  cell  operating 
temperatures  for  a  fixed  methanol  feed  concentration  of  0.5  M  at 
a  methanol  feed  flow  rate  of  1.36  ml  min-1.  In  these  experiments, 
the  polarization  curve  show  a  sharp  drop  in  the  voltage  which  is 
due  to  methanol  mass-transport  limitations.  The  current  density 
at  which  the  extrapolated  cell  voltage  tends  to  zero  is  taken  as  the 
experimentally  determined  limiting  current  density  and  is  used  for 
comparison.  It  can  be  observed  that  the  limiting  current  densities 
increase  with  increasing  methanol  feed  concentrations  (over  the 
range  of  concentrations  used  here)  and  operating  temperatures. 
Increasing  limiting  current  density  with  increasing  cell  tempera¬ 
ture  can  be  attributed  to  enhanced  methanol  mass  transport  across 
the  anode  compartment.  A  different  dependence  of  methanol  feed 
concentration  on  limiting  current  behaviour  can  be  observed  with 
high  methanol  feed  concentrations  and  this  will  be  discussed  later. 
Since  only  a  single-phase  mass  transport  model  is  used  here,  the 
model  cannot  be  used  to  predict  limiting  current  densities  accu¬ 
rately,  if  the  carbon  dioxide  produced  on  the  anode  side  is  not  in  a 
fully  dissolved  state. 

3.2.  Methanol  crossover 

A  non-zero  value  of  methanol  concentration  at  the  ACL|PEM 
interface  can  lead  to  methanol  crossover  through  the  PEM  to  the 
cathode  side  where  it  is  oxidized  without  producing  usable  protons 
and  electrons.  It  also  competes  with  the  oxygen  diffusing  through 
the  CBL.  The  mass-transport  models  developed  for  the  anode  com¬ 
partment  and  the  PEM  are  used  to  explore  the  effect  of  various 
parameters  on  this  undesirable  methanol  crossover  phenomenon. 

Fig.  4a  and  b  illustrates  the  current  density  effects  on  methanol 
crossover  for  various  methanol  feed  concentrations  ranging  from 
0.125  to  5.0  M  at  90  and  70  °C,  respectively.  Decreasing  methanol 
crossover  with  decreasing  temperature  can  be  observed.  The  effect 
of  temperature  on  the  crossover  is  related  to  the  methanol  diffu- 
sivity  through  the  PEM.  With  all  the  other  parameters  remaining 
the  same,  decreasing  the  cell  temperature  decreases  the  diffusivity 
of  methanol  and  thus  the  crossover  flux  through  the  PEM.  Fig.  4a 
and  b  also  show  that  the  dependence  of  methanol  crossover  on  the 
cell  current  density,  as  a  function  of  methanol  feed  concentration,  is 
not  unique.  At  low  feed  concentrations,  the  crossover  flux  decreases 
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Fig.  3.  Limiting  current  densities  predicted  using  the  present  mass  transport  model 
(solid  line)  and  experimental  data  (markers)  [30]  at  a  methanol  flow  rate  of 
1.36  ml  min-1  for  (a)  various  methanol  concentrations  at  a  cell  temperature  of  90  °C 
and  cathode  side  pressure  of  2  bar,  (b)  for  various  cell  temperatures  at  a  methanol 
concentration  of  0.5  M  and  at  a  cathode  side  pressure  of  2  bar,  and  (c)  for  various 
methanol  concentrations  for  a  temperature  of  70  °C  and  cathode  side  pressure  of 
3  bar. 

as  the  current  density  increases.  However,  as  the  feed  concentra¬ 
tion  increases,  the  methanol  crossover  flux  remains  almost  constant 
or  increases  with  increasing  current  density.  Thus  an  optimization 
needs  to  be  done  between  increasing  current  density  and  decreas¬ 
ing  crossover  flux. 

Methanol  crossover  across  the  PEM  takes  place  via  diffusion  and 
convection,  the  latter  arising  from  the  driving  forces  of  potential 
and  pressure  gradients.  For  low  methanol  feed  concentrations,  the 
potential  gradient  driven  electro-osmotic  methanol  crossover  flux 
is  very  small  (see  Fig.  4c),  and  as  the  current  density  is  increased, 


Fig.  4.  (a)  Predicted  methanol  crossover  flux  vs.  cell  current  density  for  various 
methanol  feed  concentrations  (0.125,  0.25,  0.5,  1.0,  2.0,  3.0,  4.0,  and  5.0  M)  at  a 
methanol  feed  flow  rate  of  1.36  ml  min-1 ,  at  90  °C,  and  at  an  air  feed  cathode  pres¬ 
sure  of  2  bar.  (b)  Predicted  methanol  crossover  flux  vs.  cell  current  density  for  various 
methanol  feed  concentrations  (0.125,  0.25,  0.5,  1.0,  2.0,  3.0,  4.0,  and  5.0  M)  at  a 
methanol  feed  flow  rate  of  1.36  ml  min-1 ,  at  70  °C,  and  at  an  air  feed  cathode  pres¬ 
sure  of  2  bar.  (c)  Predicted  electro-osmotic  methanol  crossover  flux  vs.  cell  current 
density  for  various  methanol  feed  concentrations  (0.5,  1.0,  2.0,  3.0,  4.0,  and  5.0  M) 
at  a  methanol  feed  flow  rate  of  1.36  ml  min-1,  at  90  °C,  and  at  an  air  feed  cathode 
pressure  of  2  bar. 

the  concentration  of  methanol  at  the  ACL|PEM  interface,  which  pre¬ 
dominantly  determines  the  extent  of  diffusive  methanol  crossover 
flux  across  the  PEM,  decreases  (see  Fig.  2)  and  becomes  zero  at  the 
limiting  current  values.  Therefore,  for  low  feed  concentrations  of 
methanol,  the  net  methanol  crossover  flux  monotonously  decreases 
with  increasing  current  density  and  becomes  zero  exactly  at  the  cor¬ 
responding  limiting  current  density  values.  This  behaviour  in  turn 
confirms  that,  for  low  methanol  feed  concentrations,  the  methanol 
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Fig.  5.  Predicted  diffusive,  electro-osmotic,  and  net  methanol  crossover  flux  for  a 
methanol  feed  concentration  of  0.5  M,  at  a  methanol  feed  flow  rate  of  1.36  ml  min-1 , 
at  90  °C,  and  at  an  air  feed  cathode  pressure  of  2  bar. 

crossover  flux  across  the  membrane  is  purely  diffusion-driven  (see 
Fig.  5).  However,  a  lesser  net  methanol  crossover  flux  than  the  dif¬ 
fusive  crossover  flux  can  be  observed  up  to  a  current  density  value 
of  close  to  0.1  A  cm-2.  This  is  due  to  the  negative  effect  of  pres¬ 
sure  gradient  driven  crossover  taking  place  from  the  cathode  to  the 
anode. 

For  high  methanol  feed  concentrations,  even  though  the  concen¬ 
tration  of  methanol  in  the  ACL  decreases  with  increasing  current 
density  (decreasing  diffusive  crossover  flux  can  be  seen  in  Fig.  6a), 


Fig.  6.  (a)  Predicted  diffusive  methanol  crossover  flux  vs.  cell  current  density  for 
various  methanol  feed  concentrations  (0.125, 0.25, 0.5, 1.0, 2.0, 3.0, 4.0,  and  5.0  M)  at 
a  methanol  feed  flow  rate  of  1.36  ml  min-1 ,  at  90  °C,  and  at  an  air  feed  cathode  pres¬ 
sure  of  2  bar.  (b)  Predicted  diffusive,  electro-osmotic,  and  net  methanol  crossover 
flux  for  a  methanol  feed  concentration  of  5.0  M,  at  a  methanol  feed  flow  rate  of 
1.36  ml  min-1 ,  at  90  °C,  and  at  an  air  feed  cathode  pressure  of  2  bar. 


Current  density  (A/cm2) 


Fig.  7.  (a)  Cell  polarization  behavior  predicted  by  the  present  model  (solid  lines) 
and  experimental  data  (markers)  [30]  for  various  methanol  feed  concentrations  at  a 
methanol  feed  flow  rate  of  1.36  ml  min-1 ,  at  90  °C,  and  at  an  air  feed  cathode  pressure 
of  2  bar. 

over  the  range  of  current  densities  shown,  no  significant  drop  in 
diffusive  crossover  flux  can  be  observed.  Fig.  6b  shows  that  at  a 
methanol  concentration  of  5  M,  the  net  crossover  flux  increases 
with  increasing  current  density  because  of  the  enhanced  electro- 
osmotic  crossover  flux  with  increasing  current  density  (Fig.  6b).  It 
is  expected  that  a  methanol  transport-based  limiting  current  will 
also  be  there  for  high  methanol  concentrations;  but  it  will  be  of  little 
practical  value.  The  negative  effect  of  pressure  gradient  [31  ]  is  small 
and  can  be  observed  up  to  a  current  density  value  of  0.1  A  cm-2.  For 
2  and  3  M  methanol  feed  concentrations,  the  crossover  flux  remains 
almost  constant  and  the  cell  may  reach  optimum  performance  with 
these  concentrations  with  the  set  of  operating  conditions  used 
here.  The  predicted  trends  are  similar  to  those  measured  [32]  or 
calculated  [11,33,34].  Among  the  other  parameters,  the  model  pre¬ 
dictions  show  that  the  effect  of  feed  flow  rate  on  methanol  crossover 
is  rather  slight.  However,  the  thickness  of  the  electrolyte  membrane 
has  a  significant  effect  and  this  is  in  agreement  with  the  results 
of  Jung  et  al.  [35].  These  results  show  that  the  Faradaic  efficiency 
(fuel  utilization  efficiency)  of  a  DMFC  is  considerably  reduced  when 
high  current  densities  (and  hence  high  methanol  concentrations) 
are  required. 

3.3.  Polarization  behaviour 
3.3.1.  Model  calibration 

In  order  to  predict  the  polarization  curves,  the  values  of  the 
parameters  k\,  k3  and  aa  related  to  the  anode  kinetics  (see  Eq. 


Fig.  8.  Effect  of  methanol  crossover  on  cathode  polarization  for  a  methanol  feed 
concentration  of  0.5  M  at  a  methanol  feed  flow  rate  of  1.36  ml  min-1,  at  90  °C,  and 
at  an  air  feed  cathode  pressure  of  2  bar. 
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(33))  are  needed.  Low-concentration  (0.125-0.5  M)  data  from  Ref. 
[30]  have  been  used  to  optimize  these  and  are  listed  in  Table  2a. 
It  may  be  noted  that  the  same  set  of  values  of  these  three  anode 
kinetic  parameters  have  been  used  for  all  the  results  shown  in  this 
section.  The  predicted  polarization  behaviour  for  0.125,  0.25,  and 
0.5  M  methanol  feed  concentrations  is  shown  in  Fig.  7.  It  can  be 
seen  that  a  good  agreement  (including  for  the  open-circuit  voltage) 
is  observed  with  the  set  of  transport  and  kinetic  parameters  used 
in  this  work.  Increasing  limiting  current  densities  with  increasing 
methanol  feed  concentrations  can  also  be  observed  as  mentioned 
previously.  Moreover,  it  is  worth  noting  that  the  cell  voltage  goes 
to  zero  exactly  at  the  limiting  current  densities  predicted  by  the 
mass-transport  model.  Finally,  the  effect  of  increasing  methanol 


Fig.  9.  (a)  Predicted  cathode  polarization  behavior  for  various  methanol  feed  con¬ 
centrations  (0.125, 0.25, 0.5, 1.0, 2.0, 3.0, 4.0,  and  5.0  M)  at  a  methanol  feed  flow  rate 
of  1.36  ml  min-1 ,  at  90  °C,  and  at  an  air  feed  cathode  pressure  of  2  bar.  (b)  Predicted 
crossover  overpotentials  for  various  methanol  feed  concentrations  (0.125,  0.25,  0.5, 
1.0,  2.0,  3.0,  4.0,  and  5.0  M)  at  a  methanol  feed  flow  rate  of  1.36mlmin_1,  at  90  °C, 
and  at  an  air  feed  cathode  pressure  of  2  bar.  (c)  Predicted  open  circuit  voltages  for 
various  methanol  feed  concentrations  (0.125, 0.25, 0.5, 1.0,  and  2.0  M)  at  a  methanol 
feed  flow  rate  of  1.36  ml  min-1 , 90  °C,  and  air  feed  cathode  pressure  of  2  bar. 


concentration  on  the  predicted  open-circuit  voltage  (OCV)  is  shown 
in  Fig.  7b  for  concentrations  in  the  range  of  0.125-2  M.  It  can  be 
seen  that  as  the  methanol  feed  concentration  increases,  the  OCV 
decreases,  which  is  in  agreement  with  experimental  data.  Further 
analysis  of  the  predictions  shows  that  this  is  due  to  the  increased 
cathode  overpotential  caused  by  increased  methanol  crossover  flux. 

3.3.2.  Effect  of  methanol  crossover  on  polarization  behaviour 

The  effect  of  methanol  crossover  on  the  cathode  overpoten¬ 
tial  is  illustrated  in  Fig.  8.  The  difference  between  the  polarization 
curves  with  and  without  crossover  indicates  the  voltage  loss  caused 
by  methanol  crossover  or  crossover  overpotential.  As  can  be  seen, 
methanol  crossover  results  in  serious  cathode  overpotentials  at  the 
open-circuit  conditions.  Flowever,  even  a  very  low  current  density 
can  cause  the  crossover  overpotential  to  fall  sharply  which  then 
decreases  smoothly  with  increasing  current  density. 

Fig.  9a,  b  and  c  shows,  respectively,  the  predicted  cathode  over¬ 
potential,  the  crossover  overpotential,  and  the  OCV  for  methanol 
feed  concentrations  in  the  range  of  0.125-5  M.  Both  the  overpo¬ 
tentials  increase  with  increasing  methanol  feed  concentration.  Of 
these,  the  crossover  overpotential  reduces  significantly  as  the  cur¬ 
rent  density  decreases  while  the  cathode  overpotential  remains 
high  and  increases  drastically  at  very  high  methanol  concentra¬ 
tions.  This  is  due  to  the  increased  mixed-potential  effect  caused 
by  excessive  oxygen  consumption  by  the  increased  methanol 
crossover  with  increasing  methanol  feed  concentrations.  The  pre¬ 
dicted  OCV  (Fig.  9c)  too  is  dependent  on  the  methanol  feed 
concentration  and  it  decreases  from  0.72  to  0.6  V  when  the  feed  con¬ 
centration  is  increased  from  0.125  to  2  M.  This  trend  is  consistent 
with  trends  reported  in  the  literature.  The  predicted  polariza¬ 
tion  curves  for  methanol  feed  concentrations  greater  than  0.5  M 
are  plotted  in  Fig.  10.  These  show  that  the  limiting  current  den¬ 
sity  increases  with  methanol  concentration  up  to  a  limit.  It  then 
decreases  with  increasing  methanol  concentration,  which  can  be 
attributed  to  the  mixed-potential  effect  at  the  cathode.  These  trends 
are  in  good  qualitative  agreement  with  reported  data  [36].  This  val¬ 
idates  the  present  treatment  of  the  mixed-potential  effect  on  the 
cathode  side,  namely,  the  use  of  the  effective  current  density,  which 
is  the  sum  of  the  actual  current  density  and  the  equivalent  current 
density  required  for  methanol  oxidation,  to  determine  the  cathode 
overpotential  (see  Eqs.  (59)  and  (60)). 

Finally,  Fig.  11  illustrates  the  cell  voltage  vs.  current  density  char¬ 
acteristics  for  three  methanol  feed  flow  rates.  As  can  be  seen,  the  cell 
performance  is  strongly  dependent  on  the  methanol  feed  flow  rate 
only  when  the  flow  rates  are  low.  The  increase  in  limiting  current 
density  with  increasing  methanol  feed  flow  rate  can  be  attributed  to 
the  increased  methanol  mass  transport  across  the  anode  compart- 


Fig.  10.  Predicted  cell  polarization  behavior  predicted  for  higher  methanol  feed  con¬ 
centrations  (1.0, 2.0, 3.0, 4.0,  and  5.0  M)  at  a  methanol  feed  flow  rate  of  1.36  ml  min-1 , 
at  90  °C,  and  at  an  air  feed  cathode  pressure  of  2  bar. 
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Fig.  11.  Predicted  polarization  behavior  for  various  methanol  feed  flow  rates  (1.36, 
5.36,  and  12.36  ml  min-1 )  for  a  methanol  feed  concentration  of  1.36  ml  min-1,  at 
90  °C,  and  at  an  air  feed  cathode  pressure  of  2  bar. 

ment.  The  effect  of  methanol  feed  flow  rate  on  methanol  crossover 
is  not  very  significant  for  high  methanol  feed  flow  rates,  and  the 
improvements  diminish  at  higher  flow  rates.  These  trends  are  con¬ 
sistent  with  those  reported  in  the  literature  [37]. 

4.  Conclusions 

A  comprehensive  one-dimensional,  single-phase  mathemati¬ 
cal  model  for  a  liquid-feed  polymer  electrolyte  membrane  direct 
methanol  fuel  cell  (PEM-DMFC)  has  been  developed  taking  into 
account  all  the  necessary  mass  transport  and  electrochemical 
phenomena.  Diffusion  and  convective  transport  of  methanol  is  con¬ 
sidered  in  the  ABL,  ACL  and  PEM  regions,  whereas  only  diffusional 
transport  of  the  species  is  considered  on  the  cathode  side.  The 
model  also  accounts  for  water  transport  across  the  PEM.  A  multi- 
step  reaction  mechanism  is  considered  to  obtain  a  kinetic  model 
to  describe  the  electrochemical  oxidation  of  methanol  in  the  ACL. 
Tafel  type  kinetics  are  used  to  describe  the  simultaneous  oxygen 
reduction  and  electrochemical  oxidation  of  methanol  in  the  CCL. 
The  model  fully  accounts  for  the  mixed-potential  effect  caused  by 
methanol  crossover  at  the  cathode. 

The  model  predicts  well  the  limiting  current  behavior 
(attributed  to  limiting  methanol  diffusion  on  the  anode  side)  at 
low  different  low  methanol  feed  concentrations  at  different  oper¬ 
ating  temperatures.  The  model  is  calibrated  against  experimental 
polarization  curves  and  good  agreement  is  obtained  with  other 
experimental  data.  It  shows  that  at  high  methanol  feed  concen¬ 
trations,  oxygen  depletion  on  the  cathode  side  due  to  the  excessive 
methanol  crossover  results  in  the  mass-transport  limitations.  Thus 
the  loss  of  efficiency  of  the  DMFC  due  to  fuel  crossover  at  high 
methanol  concentrations  is  captured  well  in  the  model.  Further 
model  improvements  including  multi-phase  flow,  electronic  and 
protonic  potentials,  and  water  and  heat  management  are  worth 


future  investigation  to  improve  the  predictive  capabilities  of  the 
present  model. 
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